Gravitational Waves from Braneworld Inflation 
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We investigate the generation of primordial gravitational waves from inflation in braneworld cos- 
mologies with extra dimensions. Advantage of using primordial gravitational waves to probe extra 
dimensions is that their theory depends only on the geometry, not on the microscopic models of 
inflation and stabilization. D{D — 3)/2 degrees of freedom of the free bulk gravitons are projected 
onto the 3d brane as tensor, vector and scalar modes. We found the following no-go results for a 
generic geometry of a five (or D) dimensional warped metric with four dimensional de Sitter (in- 
flationary) slices and two (or one) edge of the world branes: Massive KK graviton modes are not 
generated from inflation (with the Hubble parameter H) due to the gap in the KK spectrum; the 
universal lower bound on the gap is 1^/3/2 H. Massless scalar and vector projections of the bulk 
gravitons are absent, unlike in geometries with KK compactification. 

A massless 4d tensor mode is generated from inflation with the amplitude H/Mp, where Mp is 
the effective Planck mass during inflation, derived from the D dimensional fundamental mass Ms 
and the volume of the inner dimensions. However, Mp for a curved dS braneworld may differ from 
that of the flat brane at low energies, either due to the ilZ-dependence of the inner space volume or 
variations in the brane separation before stabilization. Thus the amplitude of gravitational waves 
from inflation in braneworld cosmology may be different from that predicted by inflation in 4d 
theory. 



PACS numbers: 04.50. -f-h, 98.80.Cq 



I. INTRODUCTION 

Higher dimensional formulation of superstring/M the- 
ory, supergravity and phenomenological models of the 
mass hierarchy has most obvious relevance to cosmology. 
Early universe inflation can universally amplify vacuum 
fluctuations of all light (minimally coupled to gravity) de- 
grees of freedom. Thus inflation serves as a powerful tool 
to probe the fields associated with the extra dimensions. 

In this paper we consider a class of models with extra 
dimensions related to the braneworld scenarios. One of 
the most interesting recent developments of the high en- 
ergy physics phenomenology is the picture of the brane 
worlds, where our 3-1-1 dimensional spacetime is the 3d 
brane embedded in the higher dimensional space. In 
application to the very early universe this leads to the 
brane world cosmology, where our 3-1-1 dimensional uni- 
verse is the 3d curved brane embedded in the bulk. A 
lot of efforts have been devoted to the search for the 
new features of the background cosmological expansion, 
inflation in the brane world settings, cosmological fluctu- 
ations which could probe the brane world scenarios and 
others. Many of these problem are rather complicated. 
For instance, the problem of gravitational waves which 
might be observed on the brane requires treatment of all 
\D{D — 3) degrees of freedom of the bulk gravitons, the 
choice of their initial state, the issue of the KK modes, 
transplanckian problem and others. 

Cosmological gravitational waves are a very useful tool 
to probe the very early universe, because their generation 
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depends on the background geometry but not on the de- 
tails of the particle physics built in the model. A historic 
example is the theory of cosmological gravitational waves 
generated from the early de Sitter stage [1] which, ironi- 
cally, was suggested even before the microscopic models 
of inflation were proposed. 

Gravitons almost freely propagate through the expand- 
ing universe, and may leave their prints on the power 
spectrum and polarization of the cosmic microwave back- 
ground. With the ever increasing accuracy of the CMB 
anisotropy experiments, it is worth to enlarge the pool of 
the primordial gravitational waves theories. 

In the familiar model of inflation in the four dimen- 
sions, gravitational waves are generated from vacuum 
fluctuations during inflation. Gravitons have two polar- 
izations, and temporal part of their wave function hk{t) 
obeys an equation which is similar to the one for the 
massless scalar field minimally coupled to gravity. 

Late time asymptotics of hk{t), which is relevant to 
present cosmological observations, contains long wave- 
length modes k. They have (almost) scale free spectrum 
k'^\hk\'^ Ri const. The dimensionless almost constant am- 



plitude 



(1) 



is defined by the ratio of the only two scales in the prob- 
lem, which are related to the gravitational waves: the 
Hubble parameter during inflation H and the Planck 
mass Mp at low energies. Usually the Hubble parameter 
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slightly varies during inflation, H = H{t), which brings 
slight tilt to the spectrum of gravitational waves (in the 
r.h.s. of (1) H H{k)). In this paper we will assume 
H = const, but will devote some discussion to the case 
when H slightly varies. 

In this paper we will study the problem of the primor- 
dial gravitational waves generation from the early uni- 
verse in the theories with large extra dimensions. As 
opposed to excitations of some other fields, gravitational 
waves propagate in the bulk. What information about 
the bulk/brane geometry can be, in principle, observable 
on the 3d branc through the gravitational excitations? 
The problem has the flavor of the well known mathemat- 
ical quest of "how one can hear the shape of a drum?" 
[2]. 

To make geometrical setting definite, we assume the 
brane world inflation. We will consider geometry of flve 
dimensional warped metric with four dimensional de Sit- 
ter slices 

ds^ =dy^ +A^{y)dsl . (2) 

The warp factor A{y) is determined self-consistently by 
the flve dimensional Einstein equations, supplemented by 
the branc boundary conditions. This formalism is briefly 
outlined in Section II. For the sake of generality, we also 
present results for a more general case of D dimensional 
warped metric with D — 1 dimensional dc Sitter slices. 
The class of metrics (2) covers many important branc 
world scenarios [3] including the Horawa-Witten theory, 
HW with with curved branes, the Randall-Sundrum (RS) 
models, RSI model with phenomenological stabilization 
of two branes, scalar sector of the supergravity realiza- 
tion of the RS model and others. The presence of a bulk 
scalar fleld (j) (a dilaton) with the potential V{(f)) or sev- 
eral scalars is often assumed. 

The warped metric (2) can be rewritten in 
"conformally-factorized" form 

ds^ = (dx' + dsl) , (3) 

where conformal factor and conformal flfth coordinate 
are deflned by 

In our treatment of gravitational waves, beginning in Sec- 
tion III, we will extensively use conformal mapping of the 
whole geometry (3), bulk and branes on its edges, into a 
simpler manifold with boundaries, with the metric 

df = n^d8'^ = dx^ + dsl. (5) 

There are several important papers on the formalism 
of metric fluctuations [4] and properties of the gravita- 
tional waves in the brane world scenarios [5-7], and will 
try to make connection with them in the course of the pa- 
per. We will work with the metric in the form (3), which 



allows separation of variables in gravitational wave equa- 
tion, as it will be shown in Section IV. Conformal map- 
ping also significantly simplifies the issue of the boundary 
conditions, as it will be explained in Section V. Just as it 
was in the four dimensional case, in D dimensional the- 
ory the graviton eigenfunctions hk,m{x) obey an equation 
for the massless scalar minimally coupled to gravity. 

In Section VI, we show that the spectrum of nonzero 
KK modes has a mass gap with the universal lower bound 

H, independently on the bulk bulk scalar fields po- 
tential V{(j)) and the presence or absence of the second 
brane. This result generalizes the finding of the | H gap 
in the Randall-Sundrum II model with a single brane in 
the AdS bulk [6]. 

In Section VII, we show that the homogeneous zero KK 
mode m = is degenerate, and can be described by only 
two physical polarization degrees of freedom, while other 
throe can be gauged by zero by the residual coordinate 
transformation. Thus, the issue of the gravi-scalar and 
gravi-vector modes in brane cosmology is resolved: they 
do not contribute. 

All together this allows us to conclude that the late 
time asymptotic of the long-wavelength gravitational 
waves generated from the branc inflation is deflned by 
the two tensor polarization of the homogeneous KK mode 
m = 0, similar to what takes place in a case of inflation 
in four dimensions. However, as we will discuss in Sec- 
tion VIII, it does not mean that gravitational waves from 
the brane inflation in all respects are identical to those 
in the four dimensional theory. Essential difference arises 
in their amplitudes. Indeed, in five dimensional theories 
there is additional scale £, associated with the fifth di- 
mension (brane separation in the HW theory, the AdS 
radius in the RS models, etc.) With three scales in the 
problem, £, H and Mp, one may have a dimensionlcss 
amplitude different from (1). Indeed, for the braneworld 
inflation, one flnds k^^^hk = j^F{Hi), with some func- 
tion F. Remarkably, the amplitude can be expressed 
economically in terms of the effective Planck mass for 
curved (de Sitter) braneworld 

k'/'h, = . (6) 

Effective Planck mass Mp^s during inflation in the 
braneworld scenario may differ from the current value 
Mp ~ 10^^ GeV, therefore prediction for the gravita- 
tional wave amplitude from brane inflation may differ 
from that from inflation in 4d theory. 

In Section IX, we will discuss implications of our results 
to the observational test of the ratio of tensor to scalar 
modes, transplanckian problem in the context of models 
with extra dimensions, and other related issues. 

We evoke the "conformally-factorized" frame (3) 
mainly for treatment of the gravitational waves. How- 
ever, conformal mapping of the brane worlds has merits 
of its own. In particular, we show that the warp geom- 
etry with branes at the edges is conformally equivalent 
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to the portion of the spacctimc with regular end of the 
world edges. For the case of flat 4d slices of (3), the world 
edges after conformal mapping are empty. 



II. BACKGROUND AND GRAVITATIONAL 
PERTURBATIONS 

In this section we outline the general formalism for 
treating background geometry with fields (such as dila- 
ton) in the bulk and on the branes, and free perturbative 

gravitational waves in this background. In general, there 
may also be scalar and vector bulk perturbations as well 
as the perturbation in the brane position. We will not 
consider them. However, we will address the question 
if bulk gravitational perturbations generate scalar and 
vector perturbations at the brane. 

The brane worlds scenario with Einstein gravity and 
scalar fields is described by the action 



S = 



1 



j y/^d^X (ii + 2£bulk)- 



(7) 




(2[/C]+2£bra„e) 



where the first term describes the bulk and the second 
term is related to the brane(s). The extrinsic curvature of 
the brane /C^,^ is defined through the normal unit vector 
riA and the tangent vierbine 



e(;.)ef^)V^"B, 



(8) 



and K. is its trace. We use "mostly positive" signa- 
ture and curvature conventions of Misner, Thornc and 
Wheeler, and denote bulk indices by capital Latin let- 
ters A,B,..., brane indices by small Greek letters iJ,,iy, 
and brane hypersurfaces by Sa, a = 1,2. Wo use square 
brackets (e.g. [fC] = /C"*" — /C~) to denote jump of the 
geometrical quantity across the brane. 
Einstein equations in the bulk are 



Rab — -^Rqab = Tab, 



(9) 



where the stress-energy tensor of the matter fields Tab 
is described by the bulk Lagrangian >Cbuik- In particular 
case of the minimally coupled scalar field with potential, 
we take £buik = -^(Vf/))^ - y((/)) and £brane = U{(j)), 
and the bulk stress-energy tensor is 



Tab = c^M,B - (5(V0)2 + V{<j>)) qab- 



(10) 



For 5d warped geometry (2) with scalar field the bulk 
Einstein equations (9) can be written as 

dy{A-'dyA) = -\{dy<l^f " ^^(0) " {A-^OyAf , (ll) 

6{A-%Af = ^{dy^f - V{ct>) + 6A-'H' , (12) 



where the second equation is a constraint. Comprehen- 
sive study of warp geometry with a scalar field based on 
the qualitative analysis of the dynamical system (11)-(12) 
was given in [10]. 

The jump of extrinsic curvature across the brane sur- 
face is related to the stress-energy tensor of the mat- 
ter contained on the brane by Israel's junction conditions, 
which in the case of scalar field with brane potential U{(p) 
are 



dU 



(13) 



From here, the assumption of Z2 symmetry around brane 
S and the warp geometry (2) implies the boundary con- 
dition for the warp factor and the scalar field 



dyA 



u 



2(D-2)' 



1 dU 

2 5^' 



(14) 



In case of the RS models we have U (0) = ±A, where the 
sign is different for the positive or negative tension A = 
const on the brane. The bulk potential is just 5d negative 
cosmological constant V^((/>) = —A. In case of the HW 
model, U{(j)) = e~'^, V{(j)) = e~'^^. Our calculations are 
valid for arbitrary f/(0) and V{(j)). As we will see, the 
results of the paper can be generalized for a more general 
case of several bulk scalars. 

Free gravitational waves obey the linearized Einstein 
equations 



5Ri = 0. 



(15) 



Linearized gravitational waves are described by the met- 
ric perturbations Hab = SgAB, and the gravitational 
wave equation is 

ahAB - '2V'^V(^AhB)C + VaVbIi + 2R'='^jJlB)C = 0, 

(16) 

where /i is a trace h = Q^^Hab- In D dimensions, grav- 
itational waves have ^D{D — 3) independent physical 
degrees of freedom. In discussions of brane world per- 
turbations, gravitational waves are often described by 
two components of the tensor modes, associated with the 
brane isometry group. However, they correspond only to 
the part of the bulk gravitational waves which have 5 
degrees of freedom. In the formalisms of [6,7] based on 
the brane isometry group, one has to consider separately 
the pieces from scalar and vector modes to complete the 
description of the 5d gravitational waves. 

In this paper we will deal with all five bulk graviton 
polarizations. We will use the "conformally-factorizcd" 
form of the background metric (3), and impose the usual 
transverse-traceless gauge for gravitational waves 



0, h = 0. 



(17) 



For this metric, it is convenient to fix some of coordinate 
freedom assuming the following orthogonality condition 
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V^CthAB = 0. (18) 

Similar coordinate gauge has been used in 4d cosmology 
where Q = i}{t) (the synchronous gauge). While the two 
conditions (17) and (18) are not necessary compatible 
for a most general background, they arc compatible for 
the brane- world metrics (2), since the orthogonality re- 
quirement is with respect to a vector = d/dx that is 
covariantly constant in the conformal frame V^w^ = 0. 

In the braneworld scenarios it is convenient to use the 
TT gauge (17), but one also can work with gravitational 
waves in other gauges. It is interesting to follow the bulk 
gravitational waves in the gauge which is usually used 
for the Kaluza-Klein compactification (however, without 
symmetry in the compact y direction). In the KK 
gauge, bulk gravitons are manifestly projected into four- 
dimensional world as four dimensional TT tensor mode 
plus scalar and vector modes. Consider the massless (ho- 
mogeneous with respect to y) free scalar and vector four 
dimensional components. After substituting perturba- 
tions into the junction conditions (13) we derive that all 
light scalar and vector projection of the free bulk grav- 
itational waves vanishes. This result will be confirmed 
further by calculations in the gauge (17) which we adopt 
in the rest of the paper. 

III. GRAVITATIONAL WAVES AND 
CONFORMAL MAPPING 

There is a useful analogy between gravitation waves 
in the the braneworld metric (3), and those in the 4d 
cosmology with the metric 

ds^ =a^{T){-dT'^ +dsl) , (19) 

where dss = jijdx^dx^ is the 3d metrics of constant cur- 
vature = 0, ±1, a(r) is the scalar factor of the universe 
and r is a conformal time. It is well known [8] that the 
tensor perturbations in this metrics can be described by 
hij which can be conformally transformed and decom- 
posed as hij = y^hk{T)Qij , where Q,, arc two-component 
transverse-traceless eigenmodes of ds|, {'V^ + k'^)Qij = 0. 
Temporal part of the eigenmodes obeys an equation sim- 
ilar to that of the minimally coupled (nonconformal) 
scalar field 

hk+(^k''-^+ 2K^ hk = 0, (20) 

where dot is derivative with respect to r. 

The warped metric (3) is a higher dimensional ver- 
sion of (19) with time and one of the spatial coordinates 
exchanged. Therefore we expect that gravitational wave 
equation (16) can be solved in the form /i^j, oc hm{x)Qij.v, 
where five gravitational polarizations are described by 
the transverse-traceless eigenmodes Q^v of 4d de Sitter 
spacetime, and hm{x) obeys an equation for minimally 



coupled bulk scalar depending on conformal fifth coordi- 
nate X- Additionally, instead of the Cauchy problem for 
graviton wavefunction hk{T) in the metric (19) we have 
Neumann boundary problem for graviton wavefunction 
hm{x) ill the metric (3). 

Let us begin with conformal transformation. All geo- 
metrical quantities and operators evaluated with respect 
to the original metric qab will be hatless, while those 
evaluated with respect to the conformally transformed 
metric qab = ^^9ab will be denoted by hat. Covariant 
derivative in physical frame will be Vj while covariant 
derivative with respect to c/ab will be V. Under general 
conformal mapping, background Ricci tensor transforms 
as 

Rab = Rab + {D- 2)f2-iVAVsr2 + 

+n-\dn)gAB -{D- i)Q-^{tnfgAB, (21) 

while its trace transforms as 

R = n^R + 2{D - i)n{an) - d{d - i){wnf. (22) 

For gravitational perturbations we have 

Kab = i^^hAB- (23) 

It can be shown that the gauge choice (17) and (18) is 
conformally invariant. Prom the wave equation in the 

physical frame (16), we obtain the gravitational wave 
equation in the conformal frame 

ahAB - 2V^V(^/lB)c + 2lf^ JlB)C 

-{D - 2)n-^n'^VchAB = 0, (24) 

where no assumptions about the background spacetime 
were made yet except the gauge (17) and (18). Although 
equation in the conformal frame appears to have addi- 
tional term compared to (16), actually considerable sim- 
plification occurs due to the simpler form of the back- 
groimd metric in the conformal frame. 

So far conformal transformation of the metric has been 
kept rather general. Now we apply the result to the con- 
formal mapping (5) with CI = f^(x). The curvature of 
the metric dP is determined by the curvature of the de 
Sitter portion, i?^,^ = {D — 2)H^ g/^v Then equation (24) 
simplifies to 

afiAB - 2V'^V(a/is)c + 2(£' - 2)H^hAB 

-{D - 2)n-^Q'h'AB = 0, (25) 

where prime denotes the derivative with respect to the 
extra dimension x- 

IV. SEPARATION OF VARIABLES IN THE 
GRAVITATIONAL WAVES EQUATION 

In this section, we continue the analysis of the gravi- 
tational wave equation (25). Since from the gauge (18) 
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it follows that = 0; only h^^i, components need to be 
considered. As = d/dx is a Killing vector for the met- 
ric in the conformal frame, the dependence on the extra 
dimension separates, and the gravitational perturbation 
can be decomposed as 



(26) 



where m is an eigenvalue corresponding to the Kaluza- 
Klein modes. The transverse-traceless conditions (17) 
are inherited by the {D — l)-dimensional tensor mode 



(27) 



while no gauge conditions are placed on extra- 
dimensional mode hm- 

Upon substitution of the mode decomposition (26) into 
the equation (25), the gravitational perturbation problem 
is reduced to two decoupled differential equations. One 
of them is equation for the {D — l)-dimensional tensor 
mode 

□g(,™) - 2V^V(^Qi7; + 2(1? - 2)H'Q<p^ = m^Q^^K 

(28) 

and the other is an equation for the extra-dimensional 
scalar mode 



h';^ + (m^ - Feff) hm = 0, 



(29) 



where effective potential Ves is determined by the bulk 
geometry 



D-2 /n" D 

77 ~ 2" T]2" 



(30) 



There is some ambiguity in the choice of the separation 
constant m^. For instance, one can include the curvature 
of the 4d slices into m^, as it is usually done for the 
gravitational waves in open or closed cosmology as in 
equation (20). Our choice of is motivated by the 
following reason: it turns out, as we will soon see, that 
in (28) and (29) has the lowest eigenvalue value = 
and is always positive. 

The equation (28) is exactly the equation for linear 
massive tensor perturbation of {D — l)-dimensional sub- 
space. Thus, the problem of gravitational waves in the 
braneworld metric (3) is reduced to two decoupled equa- 
tions, one for the massive tensor polarizations Q^^^ on 
the {D — 1) spacetime of constant curvature which we 
will analyze in detail in Section VII, and one for the min- 
imally coupled scalar mode hm{x)- These equations are 
valid for a general case of an arbitrary warp factor ^{x)y 
and for flat or curved branes of constant curvature. 

The equation (29) has a convenient "Schrodinger-type" 
structure. The effective potential (30) can be represented 
in different forms which are useful for different purposes. 
Using conformal transformation of the Ricci scalar (22) 



_R 

f22 



R + 2{D-1) 



77 



the equation (29) can be written in the form 

R 



h'L + 



hm = 0, 



(31) 



(32) 



and we have an alternative expression for the effective 
potential 



D-2 



Q2- 



(33) 



Numerical factor = \%E^ corresponds to the value 
of the coupling constant for a scalar field conformally 
coupled with gravity. This is what one expects to obtain 
for the equation of the minimally coupled (nonconfor- 
mal) scalar field after conformal transformation. Extra- 
dimensional graviton mode in the physical frame corre- 
sponds to the minimally coupled massless scalar field. 

If we return to the warp factor A{x) = l/^(x) instead 
of ^{x) in (30), master equation (29) will be written as 



h'' 



D-2 



D-AA'"^ 



hm = 0. (34) 



This equation generalizes equation for the 4d cosmolog- 
ical gravitational waves (20) (with the shifted choice of 
the separation eigenvalue k'^ + 2K). 

Let us show that there is no tachyonic modes in the 
spectrum of KK modes, i.e. is always positive. The 
equation (29) with the potential in the form (30) can also 
be written in a factorized form 



-V-V+hm = m^hm, 
where the differential operators are 



(35) 



(36) 



Simple manipulations with the equation (35) using inte- 
gration / dx between two branes, lead to the formula 



,_ J\V+hm\'dx 



(37) 



provided that boundary terms satisfy {h'!^'D^hm)\iia = 0- 
This condition is met for Z2 symmetry. 



V. BOUNDARY CONDITIONS AND 
CONFORMAL MAPPING 

The mass spectrum of the Kaluza-Klein modes is deter- 
mined by a boundary value problem in the extra dimen- 
sion, and is obtained by imposing desired boundary con- 
ditions on solutions of the equation (32). In this section 
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we discuss how the background and perturbative quan- 
tities related to boundaries are mapped under conformal 
transformation. The goal is to derive the boundary con- 
dition for the mode function h„i{x)- If Z2 reflection sym- 
metry is assumed across the brane, the constant tension 
brane is mapped to the edge of the conformal space (Jab- 
Remarkably, it turns out that after conformal mapping 
the edge regular, i.e. is not marked by the jumps in the 
curvature. This conformal mapping is sketched at Fig- 
ure 1. 








factorizable | 






1 

1 




geometry 1 
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FIG. 1. Conformal transformation simplifies geometry. 
Upper panel represents warped geometry (2) with branes car- 
rying tensions on the edges. Lower panel represents confor- 
mally-factorized geometry (5) with regular edges. 

The boundary condition for the gravitational wave 
hAB in the conformal frame at the empty Z2-symmetric 
brane is simply h'^^, — 0. Therefore, for the gravitational 
wave amplitude h„n we have a simple mirror-like (Neu- 
mann) boundary conditions 



/ J3-2 \ ' 

[p. 2 h,nj 



0. 



(38) 



In the rest of this section we shall demonstrate 
that conformal mapping of the whole geometry, the 
conformally-factorized bulk metric (3) and its brane 
boundaries, brings it to the bulk metric (5) with empty 
boundaries. This brane- world conformal mapping prop- 
erty works for a broad class of warp geometries with the 
bulk scalars with flat or de Sitter branes. 

Let us consider the effect of the conformal transforma- 
tion dp = ds^ on the surface S of co-dimension one 
embedded in the spacetime manifold, and, in particular, 
on the extrinsic curvature /C^^. Suppose the surface to be 
parametrized as = x^{(,^), then the holonomic basis 
vectors tangent to the surface E have the same compo- 
nents in both physical and conformal frames 



(39) 



Unit normal to the surface S must be scaled, however, 
to preserve the normalization 



fln^, riA = ^fiA- 



(40) 



The extrinsic curvature (8) of the surface S under con- 
formal transformation behaves as 

/c^, = n-'JCf,, ~ n-^{vn ■ h) (4i) 

with its trace transforming as 

IC^nK.-iD-l){Vn-h). (42) 

While tangential derivatives of the conformal factor 
must always be continuous across the surface for induced 
metric to match, the normal derivative need not be so. 
In particular, if the brane- world has the Z2 symmetry 
across the brane, IC~^, — —IC'^^, and 1C~^ — — and so 
the normal derivative of the conformal factor flips sign 
(VJ7-n)~ = — (VrJ-fi)^. In geometrical terms, this means 
that the gradient of the conformal factor Vfi experiences 
mirror reflection on the brane. 

Now let us consider what happens to junction condi- 
tions (13) under conformal transformation. The left hand 
side transforms as 

[Ki^ ~ ]C5ii] = niti^ - /C<5^;] + {D- 2)[Vr! • fi] 61^. (43) 

In the physical frame we have [JC^ — /C(5^] — U{(j))S^, 
where U{(f>) is the scalar field potential at the brane. On 
the other hand, the last term in equation (43) gives us 

(D - 2)[Vr! • n] = 2{D - 2)d^n\^ = -2{D - 2)%i|s, 
where A is the warp factor. As can be seen from the 
junction condition in the form (14), these two terms in 
(43) are identical, and cancel. Therefore for the junction 
condition in the conformal frame we simply have 



(44) 



This means that in the conformal frame the brane be- 
comes a regular surface, which means no special junction 
conditions or boundary terms are needed. The flat branes 
with 77 = in the conformal frame becomes an empty 
regular surfaces. 



VI. SPECTRUM OF KK MODES 

Equipped with the extra-dimensional mode equation 
(29) and its boundary conditions (38), we can analyze its 
solution ft.m(x); in particular, the spectrum of KK modes 
for the brane world warped geometries (3). 

First we notice the presence of zero mode ni^ — ho- 
mogeneous KK solution, which immediately follows from 
the equation in the form (35) 
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= 0. 



(45) 



For non-zero KK modes we have to study the eigen- 
values of the Schrodinger-type equation (33), which de- 
pend on the the shape of the effective potential Vcsix)- 
Figure 2 illustrates this function for two examples of 
the braneworld models. Left panel corresponds to the 
Randall-Sundrum models (with one or two branes) with 
negative cosmological constant A in 5D AdS bulk. Right 
panel shows VeS in the warp geometry model with bulk 
scalar field potential F(0) = Vbe~^"^, corresponding to 
the HW model. 




FIG. 2. Schrodinger potentials (red solid line) and Vcff 
(blue dotted line) for the RS model (left) and HW model 
(right) in the periodic brane-worlds between two branes. Up- 
per horizontal (green) line correspond to the gap f-ff in the 
RS model; lower horizontal (green) line correspond to the 



generic gap 



H. 



At first glance, the spectrum of massive KK modes 
crucially depends on Vcs{x)i and consequently, on the 
braneworld theory. However, it turns out that in general 
case one encounters the gap in the spectrum of KK modes 

always greater or equal to ^J^H. Exact numerical factor 

in the front of H can be slightly different (but greater) 
in different theories. 

First, let us inspect the equation (33). If the bulk 
curvature is negative, i? < 0, then the effective potential 
in the equation for the KK amplitude hm{x) has a lower 
bound 



V+h^. While h. 



obey the equation in the form (35), 
the eigenfunctions hm obey the adjoint equation 



-V- 



(47) 



In the terminology of SUSY quantum mechanics, the sys- 
tem hm is a superpartner of hm ■ The essential feature of 
adjoint representation for us is that the set of eigenvalues 
is the same for both systems hm and hm- 
Now consider equation (47). This equation again can 
be rewritten in the form of the Schrodinger equation 



VcS I h,r,. ^ 0, 



(48) 



where, using the Einstein equations (11,12) for the bulk 
scalar field, the effective potential T4ff for hm can be writ- 
ten in the form 



£1-2 , 



{D-A){D-2) 



[dyAf 



(49) 



In case of several bulk scalars (pi gradient term {dy4>)'^ in 
(49) is replaced by positive term '^i{dy(pi)^ ■ This "su- 
persymmetric" effective potential V^s for two examples, 
the RS and HW models, is plotted in Figure 2 alongside 
with VcS- From (49) one can immediately see that for an 
arbitrary scalar field theory 



> H^. 



(50) 



Therefore, for D = 5 the spectrum of KK modes has the 
mass gap of at least \I\h. 

The presence of the gap means that there will be no 
KK gravitational modes (except m = mode) which 
will contribute to the late time cosmological gravitational 
waves. This is because the allowed modes are heavy, 

TO > \J\h for _D = 5, relative to the Hubble parameter 
during inflation. As we will argue in the Section VII, 
the heavy KK modes do not survive the late-time long 
wavelength cosmologically relevant asymptotics. 



D-2 



, i? < 0, 



(46) 



and the mass gap to > H appears in the spectrum. 
This takes place in both RS models (with one or two 
branes) with the AdS background supported by the neg- 
ative bulk cosmological constant. The presence of the 
gap f-ff in the RS model with a single brane was noticed 
earlier in [6] and in the domain wall in AdS model in [5] . 

If the bulk curvature is positive, the potential Vcs can 
be smaller than (-^^2^) or even become negative. To 
analyze the mass spectrum in this, more general, case 
we will use the following trick. For a given system of 
eigenfunctions hm, let us introduce adjoint system hm = 



VII. FOUR-DIMENSIONAL MASSIVE TENSOR 
MODES 

The gravitational waves from the brane-world infla- 
tion can be factorized into the Kaluza-Klein modes hm (x) 
which obey the equation (29), and (D — 1) dimensional 

tensor modes Q^jZ\x'') which obey the equation (28). In 
this section we consider {D — 1) tensor modes. For cos- 
mological applications we will concentrate on I? = 5, i.e. 
on four-dimensional tensor modes. 

De Sitter manifold can be covered by different coordi- 
nates. It is customary in the inflationary cosmology to 
use a spatially flat dS slicing 
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= a^(r) (-dr^ + dxidx') 



(51) 



where i = 1,2,3, a(T) = —1/Ht. These coordinates 
are very convenient to study scalar fluctuations or four- 
dimensional massless gravitons. In paper [6], they were 
used to study massive transverse traceless tensor modes 
of the isometry group R? of a flat three dimensional slices 
of (51). As we already noted, this approach is not fully 
consistent, since these modes have only two polarizations 
and do not describe the five component bulk gravitons. 

To deal with the massive transverse traceless ten- 
sor modes of dS space with the bigger isometry group 
50(4, 1) we shall use the global de Sitter map 



ds^s = -dt^ + ^9 cosh^ {Ht) dnl, 



(52) 



= {-dT^ + d^l); 

COS^T ^ ' 

rfOg = + sin^ '4}{d6'^ + sin^ Odcj?). 



Equation (28) for the metric (52) is the equation for ten- 
sor eigenfunctions of the Laplacc-Bcltrami operator on dc 
Sitter spacetime. The iteration scheme to construct ten- 
sor eigenfunctions from scalar hyper-spherical harmonics 
was developed in [9]. 

First, we will consider the massless mode w? = 0. This 
is very important special case corresponding to the ho- 
mogeneous KK mode. In this case the equation (28) 
corresponds exactly to the equation for the usual four 
dimensional massless tensor mode. Tensor mode Q^^^, 
has not five but only two transverse traceless degrees of 
freedom, corresponding to the two polarizations of four- 
dimensional graviton. Technical reason for this degen- 
eracy is that there is a residual gauge freedom of the 
synchronous gauge in four dimensions which can be used 
to zero three components of tensor perturbations. When 
is non-zero, the degeneracy is lifted. 

In rn? = case solution of the equation (28) can be 
further factorized and written as [9,11] 

QS. = , Ql = -^Jnimir) T^''"(V, 0, </.) , (53) 
a(r) •' 

1 

(54) 



^2n{n+l){n + 2) 
X ((n + 2) COST - e"'^) e-'("+i)^ 



where N is normalization constant independent on 

{n,l,m), Tj"'™(?/;, 6*, (/)) arc transverse traceless tensor 
harmonics of three-sphere dU^ explicitly derived in [9], 
which have two independent degrees of freedom. 

The key point for our discussion is that in the late time 
limit t 00, T ^ 7r/2, the amplitude fnim{T) freezes out 
at the constant value fnimW) — » , ^ Similar 

treatment of massless gravitons in the planar coordinates 
(52) gives us frozen amplitude {2k)^/^ where k is 3D 



comoving momentum. Therefore KK zero modes arc, in 
principle, cosmologically observable, depending on their 
amplitude N. 

Let us turn to the massive KK modes with > 0. 
In this case there is no degeneracy of bulk graviton de- 
grees of freedom, and simple factorization (53) is not ap- 
plicable. Massive tensor KK modes in five dimensions 
are described by all five components. As we shown in 
the previous section, in the brane inflation scenario mas- 
sive gravitons are heavy, with the mass value which ex- 
ceeds the universal threshold \J^H. The key point about 

heavy gravitons is that they are not generated during in- 
flation. 

This is similar to the properties of massive scalar 
field fluctuations during inflation. Indeed, familiar scalar 
eigenfunctions of the Laplace-Beltrami operator in the 
de Sitter geometry in the planar coordinates are given in 
terms of the Hankel functions 



with the index 



{kr) 



A 



(55) 



(56) 



The physical wavelength of the quantum fluctuations 
is stretched out with time as Cosmologically inter- 

esting fluctuations have very small comoving momenta 
k ^ 0. The amplitude of interest corresponds to the late- 
time asymptotics of conformal time r ^ (or t — > oo). 
The late-time long-wavelengths asymptotic fcr — > of 
•/"fe™^ (r) crucially depends on the value of parameter 



For massless or very light modes, ^ 
tions are frozen out at the level 



H 



x/2fc3/2 



m ■ 

<C 1, the fluctua- 



(57) 



as far as its physical wavelength exceeds the Hubble ra- 



dius 



On the contrary, the heavy modes 



> do not freeze out but oscillate with decreasing 
amplitude 



rim), \ 



(58) 



as T ^ 0. 

Therefore we conclude that the heavy KK modes do 
not contribute to the cosmologically interesting late-time 
long-wavelengths asymptotics of the gravitational waves 
and scalar fluctuations. This conclusion is similar to the 
result for the scalar field fluctuations in the model with 
inflation in outer space in theories with extra dimensions 
(anisotropic inflation) with the KK compactification [12]. 
There it was shown that heavy KK modes do not con- 
tribute to cosmological fiuctuations generated from infia- 
tion. 
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VIII. BRANEWORLD DIMENSIONAL 
REDUCTION AND NORMALIZATION OF THE 
GRAVITATIONAL WAVE 

Einstein gravitational action in D dimensions is 

Sd = M^-^ J R[gAB] d'^x (59) 

plus surface terms, where Ms is D-dimensional funda- 
mental mass. If we can factorize the metric qab into 
four dimensional outer space and {D — 4) dimensional 
inner space of volume Mnnen then the action (59) will be 
reduced to the effective four-dimensional Einstein action 



^4 



^P.cff 



(60) 



where the effective four-dimensional Planck mass Mp^eS 
is related to the fundamental mass Ms via formula 



(61) 



Let us follow how this procedure works for particu- 
lar example of conformally-factorizable 4-1-1 metric (2) 
of the braneworld scenario. First we perform conformal 
mapping (22) and assume R does not depends on the 
fifth dimension x- By doing this, we restrict ourselves 
to zero-mode (with respect to x) gravitational perturba- 
tions, which is the subject of cosmological interest. By 
integrating extra-dimensional dependence out, we have 



Sd = M 



D-2 



J n^-''{x)dx J R[9ab] 



(62) 



where we have omitted derivative of Q terms for brevity. 
The effective action for the four- dimensional theory is 
then 



S4 = Ml^^s j R[g^u\yf^d^x + ..., 
where effective four-dimensional Planck mass is 



Mil 



S*off 



(63) 



(64) 



and £eff is the effective size (volume) of the inner dimen- 
sion. For braneworld scenario with compact inner dimen- 
sion X) it is 



X2 



J/2 



^-\x)dx = / A\y)dy. 



(65) 



The warp factor A{y), in principle, depends on the brane 
curvature (defined by H) through the Einstein equation 
(12). In particular, it is different for de Sitter brane ge- 
ometry with the Hubble parameter H and for the flat 
brane with H = Q. Independently, physical separation 



between curved brancs may not be the same as the sepa- 
ration of flat brancs at low energies. Indeed, in the brane- 
world cosmology energy-momentum tensor at the brane 
is time dependent, and so is the extrinsic curvature, i.e. 
the brane embedding. Therefore the effective volume ^eff 
of the extra dimension is changing in the course of cosmo- 
logical evolution between inflation and present universe. 

Thus, the effective Planck mass Mp^eff should be ex- 
pected to be different during inflation and in the present 
universe. 

Let us return to the gravitational waves. In the actions 

(59), (60) we split the curvature into the background part 
plus gravitational wave perturbations. To construct the 
graviton field theory we keep perturbations up to the 
second order. Contribution to the action from the second 
order graviton perturbations reads as 



5S. 



Ml jn-^ \h,AB;ch^''-''' + . 



^d^x. (66) 



As we found in the previous section, only zero KK mode 
contributes to the interesting cosmological effects. For 
these modes integrant in (66) does not depend on the 
fifth dimension, and we obtain 



5Sa 



M^ 



P.eff 



/ 



1/, hf^^\P A_ 



-gd'^x. (67) 



Thus, we reduce the problem back to the four- 
dimensional theory of gravitational waves from infla- 
tion. The standard approach is to rewrite action (67) 
in terms of scalar field ip with the canonical kinetic term, 

h, and quantize the scalar field (p in the dc 



V2 



Sitter background. Scalar field vacuum fluctuations are 
amplified by inflation to the magnitude k^^'^^pk = (see 
also (57)). As a result, we obtain that long-wavelength 
gravitational waves are generated from inflation with the 
scale free spectrum and the amplitude 



H 



M, 



P,eflF 



(68) 



Thus, the only difference with the results for the gravita- 
tional wave amplitude in the usual four-dimensional the- 
ory is the effective Planck mass Mp^eff instead the flxed 
Mp, namely 



H 



Mp\ £in{ 



(69) 



where inow, ^inf arc the effective sizes of extra dimension 
at the inflationary stage and today. 

Let us discuss how Mp^eff at the stage of inflation in 
the braneworld scenario may be different from the cur- 
rent Mp. It is instructive to consider the simple Randall- 
Sundrum model with a single inflating brane in flve di- 
mensional AdS bulk. There we have the conformal factor 
^(x) = ^TF^- Then exphcit integration in (65) gives 
[6] 
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(70) 

For the ratio involved in formula (69) we have the 
asymptotics 

^now f 1 foriJ^ < 1 
^i„f ^ \ ^{H£) for i/^ > 1 • 

If brane inflation occurred at scales smaller than the cur- 
vature of the extra dimension (i.e. H£ » 1), and today 

universe is almost flat (i.e. H£ <C 1), the gravitational 
wave amplitude generated by inflation is enhanced 

compared to our expectation from standard cosmology. 

Another example where Mp.cff is changing in the 
braneworld scenario was considered in [7]. It can be 
shown that the model of [7] effectively corresponds to the 
case of two parallel inflating branes separated by scale £ 
embedded in flat empty five dimensional space. There it 
was assumed that H£ <C 1. 

We are not aware of convincing arguments why H£ 
cannot be of order of unity or larger. The case of Hi <^ 1 
corresponds to two (almost) parallel branes separated by 
the scale ^ £ which is small compared to their curvature 
radius H~^. It looks natural to have as initial conditions 
two branes parallel within the casually connected region 
[13] . From point of view of four dimensional observer this 
is the region ~ However, casually connected region 

in five dimensional space is defined not by H but by five 
dimensional curvature. For instance, in the RSI scenario 
with the AdS bulk already for H£ ~ 1 the function £cs 
in (70) drops by factor 2 (it is a very steep function of 
H£). The case H£ > 1 gives us interesting possibility to 
have significant departure from the standard prediction 
of the inflation in the four dimensional theory. 

So far we discussed the amplitude of the gravitational 
waves at the stage of inflation. Will formula (68) be 
the flnal answer for the gravitational wave amplitude in 
the late time cosmology? The reason for this question 
is the following. The amplitude hk may somehow be 
affected during cosmological evolution between inflation 
and the latest stages because Mp^eff niay change its value 
between inflation and current value (say, after stabiliza- 
tion after inflation).* Although equation for free grav- 
itational waves does not depends on the gravitational 
coupling, situation with the varying Mp is more deli- 
cate. One may model this situation in framework of the 
Brans-Dicke theory, where BD scalar has two constant 



*We thank Shinji Mukohyama for this comment. 



but different asymptotics. Simple estimations with this 
approach indicates that there is no additional variation 
of the amplitude hk- This issue deserves further investi- 
gation. 

IX. DISCUSSION 

Among the many features of superstring/M theory, the 
existence of extra dimensions has the most direct rele- 
vance to cosmology. Gravitons propagate in all - outer 
and inner - dimensions. Bulk gravitons in D dimensions 
have ^D{D—3) components. Only two of them appear as 
a transverse traceless tensor mode from the point of view 
of a four dimensional observer, while the others are pro- 
jected into the 3-1-1 dimensional world as (gravi)scalars 
and (gravi)vectors. 

Vacuum fluctuations of all light (minimally coupled) 
degrees of freedom are unstable during inflation and ap- 
pear as long-wavelength perturbations in late time cos- 
mology. Thus, outer space inflation acts like an amplifier 
and stretcher of light modes. Using inflation as a tool, 
we may probe extra dimensions by examining light modes 
associated with extra dimensions which may be produced 
by inflation. 

It is especially appealing to investigate gravitational 
waves in the context of extra dimensions plus inflation. 
This is because the theory of small gravitational fluctu- 
ations depends on the underlying geometry but not on 
the microscopic model of its realization. Additionally, 
gravitons propagate almost freely throughout the evolu- 
tion of the universe, thus carrying information about the 
primordial universe. 

There are two distinct types of scenarios with com- 
pactification of extra dimensions from the point of view 
of the theory of primordial gravitational waves. One type 
is Kaluza-Klein compactiflcation. The theory of gravita- 
tional waves from inflation in models with KK compact- 
iflcation will be considered separately. Another type of 
model arises from the braneworld scenarios with heavy 
branes (including a single brane in AdS bulk or compact 
warped inner space in between two branes). In this paper 
we investigated primordial gravitational waves in generic 
geometries of the warped inner space with outer space 
(brane) inflation. 

We found the following no-go results for gravita- 
tional waves from inflation in the generic case of these 
braneworld scenarios: 

Massless mode of bulk gravitons is projected to the 
usual four dimensional two component TT tensor mode, 
while their scalar and vector projections are absent. This 
is different from KK compactiflcation where light bulk 
gravitons are projected to tensor, vector and scalar mass- 
less components in the four dimensional world. 

Massive KK modes are not generated from inflation 
due to the gap in the KK spectrum. The gap is present 
in all models including models with two branes. The 
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universal lower bound of the gap is 



Am > 



H. 



(73) 



This result generalizes the finding of the gap (of magni- 
tude |-ff) for the RS model with a single brane embedded 
in AdS [6]. 

Among all massless projections of the bulk gravitons, 
only the four dimensional TT massless modes are gen- 
erated from brane inflation. As in the standard four di- 
mensional theory, their spectrum is scale free, but the 
amplitude is given by the formula 



H 



P,inf 



(74) 



where Mp^i^t is the four dimensional Planck mass for the 
curved (de Sitter) observable brane. Suppose the inner 
dimensions are stabilized after inflation (say due to SUSY 
breaking). This is not too radical an assumption bear- 
ing in mind that the usual energy scale of inflation is 
H ~ 10^^ GeV and the usual energy scale of SUSY break- 
ing is ~ 10^^ GeV. Therefore theoretically it is not nec- 
essary that Mp.inf during inflation coincides with present 
day low energy Mp. The difference may come from vari- 
ation of the size of the inner dimension(s) (i.e. the limits 
of integration in formula (65), or from an imprint of the 
brane curvature H on the warp factor A{y). A change 
in Mp eff by a factor of 2 is already interesting. Thus in 
general, the connection between gravitational wave am- 
plitudes and energy scale of inflation H in theories with 
extra dimensions may be different from that in four di- 
mensional theory. 

Most important test of the primordial cosmological 
fluctuations from inflation is the ratio between the ampli- 
tude of tensor and scalar perturbations T/ S. To calculate 
scalar perturbations, one needs to specify the model of 
inflation (say in terms of the scalar field potential at the 
brane). Already the simplest model of inflation in the 
RSI scenario without bulk scalars predicts that T/S may 
be different from that predicted in the four dimensional 
theory. The theory of scalar perturbations in braneworld 
scenarios with bulk scalars is even more involved. In gen- 
eral, we expect that T/S from braneworld inflation will 
be model-dependent and depend on many parameters. 

Finally, we will make a comment on the transplanck- 
ian problem in the context of extra dimensions plus in- 
flation. The idea beyond this problem is the following 
[14]: In an expanding universe the physical momentum 
of quantum fluctuations is redshifted to smaller values. 
Going backwards in time, we have to deal with ever in- 
creasing momenta which at some instant will be equal to 
the fundamental mass scale Ms- At this point and be- 
yond the quantum field theory approach to fiuctuations 
in an expanding universe (say during inflation) should be 
replaced by a string theory description, which is not yet 
available. For quantum fluctuations from inflation in the 



4d theory there is, however, a convenient phenomenolog- 
ical coding of the problem in terms of the Bogolyubov 
coefficients ak,f3k for the initial vacuum state: instead of 
the vacuum eigenmode /fe(r) at inflnity (given in terms of 
the Hankel function (55)) one can use the mode function 



fk{T) ^ afcA(T) +/3fe.4*(r), \ak\^ - \Pk\ 



(75) 



and discuss further the effects from (3^ , if any. If the sizes 
of the extra dimensions are larger than Mg^ which takes 
place in many phenomenological schemes (like the HW 
theory), then we must simultaneously deal with extra di- 
mensions and the transplanckian problem. It is, however, 
straightforward to combine the formalism of cosmologi- 
cal fluctuations and that of the transplanckian problem. 
If the sizes of the extra dimensions are large enough to 
work there in the supergravity approximation, we can use 
the quantum field theory of fluctuations in the bulk, as 
we did in this paper. Then the transplanckian story will 
again be coded in terms of (3k for bulk quantum eigen- 
modes (of gravitons, for instance). 
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